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We analyze the non-Markovian dynamics of a quantum system subject to spontaneous 
collapse in space. After having proved, under suitable conditions, the separation of the 
center-of-mass and relative motions, we focus our analysis on the time evolution of the 
center of mass of an isolated system (free particle case) . We compute the explicit expression 
of the Green's function, for a generic Gaussian noise, and analyze in detail the case of an 
exponential correlation function. We study the time evolution of average quantities, such as 
the mean position, momentum and energy. We eventually specialize to the case of Gaussian 
wave functions, and prove that all basic facts about collapse models, which are known to be 
true in the white noise case, hold also in the more general case of non-Markovian dynamics. 
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INTRODUCTION 



Among the several models of spontaneous wave function collapse which has been considered 
so far, the so-called QMUPL (Quantum Mechanics with Universal Position Localization) model 
is particularly interesting, being it a very good compromise between mathematical sirnplicity and 
physical adequacy. It was first introduced by Diosi [1|, |2] and subsequently studied in 
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sirnplicity and 



both from the mathematical as well as physical point of view. It is particularly relevant 
because it is the simplest model describing the evolution of the wave function of a system of N 
distinguishable particles, subject to a spontaneous collapse in space; as such, it can be analyzed in 
great mathematical detail. The model is defined by the following stochastic differential equation 
(SDE): 

N ^ N 

■^Hdt + '^ \/Ki{qn - {qn)t)dWn,t ~ 2 ^ ^niqn - {qn)tfdt 



dM{^}) 



n=l 



n=l 



M{^}); (1) 
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where the symbol {x} = xi, X2, ■ ■ ■ Xn denotes the coordinates of the particles (for simplicity, we 
will work in one spatial dimension). The operator H is the standard quantum hamiltonian of the 
composite system; g„ is the position operator associated to the n-th particle and {qn)t = (V'tknIV't) 
denotes the quantum expectation value of Qn- The stochastic processes Wn,t are independent 
standard Wiener processes defined on a probability space {0,,J^,F), and the parameters A„ are A'^ 
positive coupling constants which is convenient to take proportional to the mass of the particle 
according to the formula |9|: 

An = — - Ao, (2) 

mo 

where rUn is the mass of the n-th particle, mo is a reference mass which, at the non relativistic 
level, is reasonable to take equal to the mass of a nucleon (mo — 1.67 x 10^^'' Kg), while Aq is the 
only true new parameter of the model, whose value sets the strength of the collapse mechanism. 

The numerical value of Ao has to be chosen in such a way that: i) the model reproduces 
quantum mechanical predictions for microscopic systems; ii) it rapidly induces the collapse of the 
wave function describing the center of mass of a macroscopic object. In the literature, two quite 
different values for Ao have been proposed, the first by Ghirardi, Rimini, and Weber [12], and the 



second by Adler 



A^^"^ ~ 5.00 X 10-3 m-2 sec-\ (3) 
^Adicr ^ ^ -^2 X 10^=^2 m"2 ggp-i^ (4^ 

GRW's choice is motivated by the requirement that superpositions of order 6.02 x 10^^ nucleons, 
displaced by a distance of at least rc = cm, be suppressed within 10"'^ sec. Adler's choice 
instead is motivated by the requirement that the collapse occurs already at the level of latent 
image formation. More specifically, GRW set Agrw — 10~^^ sec"^, where A prw is the collapse rate 



of the GRW model |l^; Adler instead set 7csl ^ 2 x 10 ^^^^ cm^ sec ^ [ij], where 7csl is the 
noise-strength coupling constant of the CSL model; the relation between Ao of our model and Agrw 
and 7csL is BQ: Ao = aAGRw/2 = a^/^-fcs^/lGny^ with a = l/r^. 

Eq. ([1]) is manifestly not linear, which makes it difficult to analyze, in particular it makes it hard 
to find its solutions. The way to circumvent this obstacle is to linearize the equation, according to 

n 

the following prescription [S]. Gonsider the linear SDE: 



dM{x}) 



^ N ^ N ' 



2 

n=l n=l 



Mix}), (5) 
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where the stochastic processes ^n,t are N independent standard Wiener processes with respect to a 
new measure Q. It can be shown that HV'tlP is a martingale, which can be used to generate a new 
measure from a given one: the measure Q introduced here above is chosen in such a way that |P 
is the Radon-Nikodyn derivative of P with respect to Q, i.e.: dF = H^tpfiQ. Moreover, Girsanov's 
theorem |l5| states that the two sets of Wiener processes {Wn,t} and {S,n,t} are related as follows: 
dWn,t = din,t - 2^/\{qn)tdt. 

Given these ingredients, it is easy to relate the solutions of Eq. ([5]) to those of Eq. ([1]): given a 
solution (j)t of Eq. ([5]), one first considers the normalized state V't = </'t/||0t ||) and then replaces the 
noises {S,n,t} with {Wn,t} according to the formula given above. It is not difficult to show that the 
wave function so obtained solves Eq. ([1]). Further details can be found in j^. 

The aim of this paper is to analyze the generalization of the QMUPL model to non-Markovian 
Gaussian random processes. As discussed in , [3] ) the generalization of the linear Eq. ([5]) to 

the non-Markovian case takes the form: 

N N ,t 



-H + y J XnqnWnit) - 2 } 

n _ /n OWm(S 



mx]). (6) 



n=l n,m=l 

where the noises Wn{t) are now supposed to be Gaussian random processes whose mean and cor- 
relation function, expressed with respect to the measure Q, are equal to: 

EqlWnit)] = 0, EQ[Wnit)Wm{s)] = Dnm.{t,s). (7) 



The correlation function is assumed to be real, symmetric and positive-semidefinite [li 

Like in the white-noise case, Eq. ([6]) does not preserve the norm of the wave function. Accord- 
ingly, we assume that the physical states are the normalized states: 

m^) - (8) 

moreover, we assume that the physical probability measure is P, which is related to the measure 
Q, by means of which the statistical properties of the noises Wn{t) have been defined, according to 
the formula: 

dF = UtfdQ. (9) 

Measurable quantities are given by expressions of the form: Ep[(^/^f|0|'(/'t)], where O is a suitable 
self adjoint operator. According to relation ([9]), the following equality holds true: 

Mii'tmt)] = Ep[(</.f|o|<At)]. (10) 
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We then have the following useful results: physical properties can be computed directly from 
the solution of the linear equation Q, if one averages with respect to the measure Q, instead 
of the physical measure P. For this reason, in the following we will focus our attention only on 
Eq. (I6|) , without trying to investigate the form of the corresponding non-linear but norm-preserving 
equation. This will be the subject of future analysis. 

The goal of the paper is to unfold the dynamics described by Eq. ([6]). We will focus our attention 
on the free particle case, which can be solved exactly. The results can be generalized to the case of 
harmonic oscillators, and in general to any equation of the type containing expressions which 
are at most quadratic in the operators q and p. For other type of systems, perturbation techniques 
can be employed 



201, 
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221. We list the main results we have obtained. 



Result 1 (Sec. H]). Under the conditions: 



N 

E 

m=l 



DmN{t, S) 



0, 



n 



l,2,...iV-l, 



(11) 



the center-of-mass motion and the relative motion decouple. In particular, when Dnm{t,s) = 
^nmD{t, s), then the two equations for the center-of-mass wave function c/)™' and relative-motion 



wave function 



are: 



dt 



+ VAQ?i;(i) - 2VaQ / dsD{t,, 
" Jo 



6w{s) 



(12) 



di 



N 



N-1 



.nm I dsD{t,s) 

/o 6Wm{s) 



if, 



(13) 



n=l n,m=l 

where ■w{t) and A are defined in ([32]) and ([39|) . and idn{t), A„ are defined in ([^5]) and (06]). 
Result 2 (Sec. IIIip . The Green's function associated to a free particle of mass m (H = p^/2m) is: 



G{x,t; xo,0) 



m 



■ exp 



a/a /■* 

2^ (xoz'(O) - xz'(t)) +^J^ dsw{s)z{s) 



(14) 



2i7rhtu{t) 

where z{t) solves Eq. (|77p . with boundary conditions: z(0) = xq, z{t) = x, while u{t) is defined in 
Eq. (I82p . In particular, under the condition D{t, s) = D{\t — s\), which represents time translation 
invariance of the noise w^t), the Green's function (|14p has the simpler structure: 



G{x,t;xo,0) 



m 



2i7rhtu{t) 



exp [-^(a^o + x"^) + BtXQX + CiXq + Vtx + £i 



where: 



At = kfm. St = 2km, 



k :-- 



ih 
2m 



(15) 



(16) 
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are deterministic coefRcients, while: 



Ct = -kh't{0) + ^l dlw{l)ft{l), (17) 

Vt = kh[{t) + ^ f dl w{l)Mt - I), (18) 
^ Jo 

£t = ^ f dlw{l)ht{l), (19) 





are random coefficients. The function /(s) is solution of Eq. ()84p . with boundary conditions 
/(O) = 1, f{t) = 0, while h{s) is defined in Eq. (jQSp . As we see, in the non-Markovian time 
translation invariant case the free-particle propagator has the same structure as in the white noise 
case. 

Result 3 (Sec. llVj) . Using the Green's function (jl4p . one can rewrite the non-Markovian SDE for 
a free quantum particle in a simpler way, replacing the functional derivative with combinations of 
the position and momentum operators. The collapsing equation takes the form: 

d 



h 2m 



+ \fXqw{t) - 2\q [ dsD{t, s){qat{s) + pbt{s) + Ct{s)) 
Jo 



Mx) , (20) 



where at{s), bt{s), q(s) are defined in (|107p - (|109p . This result agrees with the Ansatz first intro- 
duced in Sec. 4.3 ofRef. Q. 

Result 4 (Sees. IVl and rVlj) . In the white noise limit, one recovers well-known results in literature. 
More interestingly, we have solved explicitly the case of an exponential correlation function: 

D{t,s) = |e'^l*-^l (21) 

in which case the functions /(s) and h(s), through which the coefficients At~£t of the Green's 
function can be computed, are: 



J2k [''t sinhffc(t — s) + cosht;fc(t — s) 



^ [2c + sinht;fct + ti^ cosh t;fet] ' ^ ^ 



while h{s) is given by ([USp with 



h^^s) = / Ml){w"{l)-^Ml))dl, (23) 



m 



fsH) = (24) 
^ k 

Result 5 (Sec. IVII|) . We have proven the so called imaginary noise tricky which we use to compute 
the average value of physical quantities, in particular the average position q, average momentum 
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p, and average kinetic energy Hq = p /2m: 

^EQ[{q)t] = ^EQ[{p)t] , (25) 

jEQ[{p)t] = 0, (26) 

-EQ[{Ho)t\ = / dsD{t,s). (27) 

The above equation prove that a free quantum particle moves, in the average, hke a classical 

particle, but its kinetic energy is not conserved, not even in the average. This is an expected 



feature of the model 



2l|. 



Result 6 (Sec. IVIIB . We have studied the time evolution of Gaussian wave functions, analyzing 
in particular the time evolution of the spread in position, and how the fluctuations around the 
average of both the mean position and mean momentum scale with the mass of the particle. In all 
cases, we have recovered the results which hold true for the white noise models. 

Before concluding this introductory section, two comments are at order. First, it may seem 
rather remarkable that we have been able to compute the Green's function associated to a non- 
Markovian equation. In fact, dynamics with memory terms in general do now allow for solutions 
expressed in terms of a Green's function, whose nature is strictly linked to the Markovian character 
of the evolution. The fact that in our case we have been able to compute the Green's function 
appears less surprising if one looks at the way in which Eq. ([6]) has been derived: originally, the 
non-Markovian evolution was presented only in terms of a propagator (see the pioneering works 



of Feynman and Vernon on decoherence [23!] and of P. Pearle [2J, |25| on collapse models) and 
only afterwards a non-Markovian equation was devised, having that propagator as the associated 
Green's function [l^. [2^. 

Second, it may appear odd that the integration interval in Eq. Q starts from and not from 
—00, as required by a truly non-Markovian dynamics. Here we are implicitly making the assumption 
that the state at time s = suffices to compute the subsequent evolution, without any need to know 
also the previous history of the system. This assumption, which actually is an approximation if the 
equation is taken as a fundamental equation, is based on the following argument. We assume that 
before s = the system had enough time to reach some equilibrium state, which does not depend 
on the way it has been reached; this is the physical argument behind the idea that the history prior 
to s = is unimportant. Then, at time s = a sudden change in the system occurs — e.g., due to 
a measurement-like interaction: Eq. ([6]) describes how the system evolves thereafter. 
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II. SEPARATION OF THE CENTER-OF-MASS AND RELATIVE DEGREES OF 

FREEDOM 

Let us suppose that the total Hamiltonian H of Eq. ([6]) separates in a term referring only to 
the center of mass and one referring to the relative degrees of freedom: H = H^^ + H^^i. In this 
section we analyze under which conditions the dynamics of the two types of motion decouple from 
each other. To this end, let us introduce the center of mass position operator Q and the relative 
position operators qn defined as follows: 

^ N N 

Q = J-f'^^nqn: M = ^nin, (28) 

n=l n=l 

qn = qn- Q- (29) 

We will choose Q and qi,q2, ■ ■ ■ , qN-i as a set of independent variables, while (Jn is given in terms 
of the other operators according to the formula: 

N-l 

i 

qN 



y^mnqn- (30) 



rriN , 

n=l 



The second term on the right hand side of Eq. ([6]) trivially separates in a term depending only 
on Q and a term depending only on the relative variables, as the following equality holds true: 



N N N-l / X \ 

VKqnWnjt) = Q'^y^Wnit) + y%^qn \Wn{t) - ^WN{t)j ■ (31) 

n=l n=l n=l ^ ^ 

The above expression suggests to define the following new stochastic processes: 

1 ^ 

w{t) = — Y^y^wnit), (32) 

^ n=l 

Wn{t) = Wnit) - ^WN{t), n = 1, 2, . . . , iV - 1, (33) 

Aiv 

with A yet to be defined. The first noise is associated to the center-of-mass motion, the remaining 
ones are associated to the relative motion. In such a way, Eq. ()3ip becomes: 

N N-l 

VKqnWnit) = V\Qw{t) + y^ \fKiqnWn{t). (34) 
71=1 n=l 

The new noises here introduced are still Gaussian processes with zero mean; a crucial role is of 
course played by their second moments. The correlation function between w{t) and Wn{t) is equal 
to: 



1 ^ 

^ m=l 



^mDmn{t, s) — DmN{t, s) 



n = 1,2,... N-l. (35) 
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As we shall see soon, the necessary condition in order for the two types of motion to decouple is 
that the above expression vanishes 



N 

E 

m=l 



■DmN{t,s) 



0, 



n = l,2,...iV- 1, 



(36) 



i.e. that w{t) be statistically independent from the other noises. This condition is automatically 
satisfied if the original noises Wn{t) are independent and identically distributed, i.e. if 



DnmitiS) — 6nmD(t,S^, 



(37) 



which is what it is usually assumed in white-noise collapse models (see Eq. ([5])). Throughout the 
rest of the section we will assume that this condition is satisfied. 
The correlation function of w{t) is equal to: 



EQ[u;(t)u;(s)] 



(38) 



this expression suggests to set the value of the coupling constant A first introduced in Eq. ([32 
equal to: 

M 



A 



n=\ 



mo 



Ao) 



(39) 



so that the correlation function of the center-of-mass noise is identical to that of the noises as- 
sociated to each single constituent. Thus, in accordance with what happens in the white noise 
case, the collapse coupling constant associated to the center-of-mass motion is proportional to the 
total mass of the system; as we shall see in the subsequent sections, this means that the collapse 
mechanism scales with the size of the system. 

According to the previous assumptions, the correlation functions of the noises Wn{t) are equal 

to: 



^Q[Wn{t) Wm{s)] = Cnm.D{t,s), 



Cn 



+ 



(40) 
(41) 



Aat 

with n, m = 1, 2, . . . , — 1. It easy to prove that the matrix C, having Cnm as coefficients, is 
positive-definite. In fact one can write C = I — XY^, where X, Y are two A^ — 1-dimensional 
vectors, Y = —X and 



X 



1 

/Ajv 



(42) 



According to Statement 3.5.14 at page 93 in [27*], C is positive definite iff X'''Y < 1, inequality 
that can be easily verified with our definitions of X and Y: 

1 A 

^"Y = -i::EA. = i-Tr<i- (43) 



1=1 

Moreover, C is real and symmetric, thus there exists an orthogonal matrix O which diagonalizes 
it: 

O^CO = D, Z) = diag((ii,...div-i), (44) 

where the eigenvalues di of the matrix D are real and positive, since C is positive definite. 

It is convenient to "diagonalize" the processes . . . , WN-i{t) by defining the new stochastic 

processes: 

1 

wn{t) = -7= Y^Ol^Wmit); (45) 



^ " m=l 

an easy calculation shows that they are independent and identically distributed Gaussian processes 
with zero mean and correlation function equal to D(t,s). By defining also: 

Qn ~ rr — ^ ^ ^nm \l^rn Qm ; '^n — dj^ , (46) 



V m=l 

one can easily verify that the second term on the right hand side of Eq. ()34p can be rewritten as: 

N-l N-l 



VK.qnWn{t) = \/\iqnWn{t) . (47) 

n=l n=l 

To summarize, the second term on the right hand side of Eq. ([6]) becomes 

N N-l 

y/K,qnWnit) = VXQw{t) + VX^qnWnjt) . (48) 

n=l n=l 

We now examine the third term of Eq. ([6]). After some easy calculations, one finds the following 
result: 



^-1 •''0 



5Wm{s) 



n,m=l 



VXQ [ dsEQ[w{t)w{s)]—^ + VXQy" [ dsEQ[w{t)wnis)] ^ , 

Jo ^UJ{S) ^ Jo dWn{s) 

N-l t ^ N-l t ^ 

\/A^9n / ds¥.Q[wn{t) w{s)] —-— + y/K.qn / lis Eq [w^ (t) (g)] (49) 
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As we see, the second and third term at the right hand side of the above equation vanish when 
condition ()36|) is satisfied, in which case the variables relative to the motion of the center of mass 
and those relative to the motion of the relative degrees of freedom decouple. This is the desired 
result. Moreover, under condition ([37|) one can rewrite Eq. (jl9]) as follows: 

N 



Jo SWm{s) 

n,m=l 

= VXQ / dsD{t,s)——-+ V ^/Xnqn 

Jo ^^(■5) .t:li Jo 



6w{s) 

VXQ [ dsD{t,s)-^^ 
Jo ^w{s) 



* 6 

nHn I ds CnmD{t, s) 



n,m=l 
N-1 



6w^ 



m \ 



+ V VA^9„ / dsD{t,s)—- 
„ Jo ^Wni 



(50) 



To summarize the results so far obtained, we have seen that when condition (j36p is satisfied, the 
center of mass motion and the relative motion decouple; under the simplifying assumption ()37p . 
the equation for the center of mass becomes: 
d 



dt 



4^c.n + VAQu;(t)-2VAQ [ dsD{t,s)-^^ 
h Jq dw{s) 



(51) 



where wit) and A have been defined in (j32p and (139p respectively. The equation for the relative 
motion instead becomes: 

N N-l 



d_ 

IE 



't 



n=l 



n,m=l 



■■nm I dsD(t,s)-^ — -r— 
lo OWm[S) 



(52) 



where Wn{t) and A„ have been defined in (j45p and (j46p respectively. 

We stress once more a crucial feature of the model, the amplification mechanism: according to 
Eq. ([39]) . the coupling between the noise and the center of mass, thus the collapse strenght, scales 
with the size of the system. On the other side, relative degrees of freedom are coupled to the noise 
through A„, which do not scale with the size of the system, and remain small. This is the way in 
which the model described by Eq. ^ describes both classical properties of macroscopic objects 
(large values of m) and quantum properties of microscopic systems (small values of m). 



III. SOLUTION OF THE EQUATION FOR A FREE QUANTUM PARTICLE 



Now we restrict our attention to the dynamics of a free quantum particle of mass m. Eq. ^ 
then becomes: 



di 



+ V\qw{t) - 2VXq / dsD{t,s)—— 
a zm Jq ow[sj 



(53) 
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where D(t,s) is the correlation function of the Gaussian noise w{t). As discussed in the previous 

section, this same equation describes also the dynamics of the center of mass of an isolated system 

(in particular, a macroscopic object) of total mass m. ^ ^ 

Aim of this section is to find the explicit solution of Eq. (j53p . In [26] it has been shown that 

the Green's function G{x, t; xq, 0) of Eq. (I53p allows for the following path-integral representation: 

rq{t)=x 
/ 

JqiO)=xo 

where, as indicated, the integration is carried out over all the paths connecting q{0) = xq to 
q{t) = X. The 'action' S[q], which contrary to the standard quantum case, has both a real and an 
imaginary part, is: 



G{x,t; xq,0) 



(54) 



ds 



'-^ q'\s) + VXq{s)w{s) - Xq{s) I drq{r)D{s, 



2h 



(55) 



That Eq. ([55]) represents the correct action associated to Eq. ([53]) can be easily verified by checking 
that the wave function 



dxoG{x,t;xo,0)(f>o{xo) , 



(56) 



solves Eq. (|53p . The advantage of the path integral representation of the Green's function, with 
respect to the standard representation associated to a differential equation like ([53p . is that it 
avoids resorting to the functional derivative of the noise, which is a source of major complications. 
We now compute the path integration in (1541) . 

n . . . . 

Following the standard Feynman polygonal approach [23], we divide the time interval [0,t] in 
N subintervals, each of length e = t/N; the intermediate time points are defined as: = ke. The 
path integral is then understood as the limit ^ oo of a multiple integral over the N — 1 variables 
Qk = q{tk),k = 1, . . . , N - 1: 



G(x, t; Xq, 0) = lim GN{x,t; xo,0) , 

Af— »oo 



with 



GAr(x,t;xo,0) 



m \ — 



2'iTihe 



N-l 



W dqk e' 



■SN[q] 



k=l 



and Sn[q] is the discretized form of the 'action', which reads: 



N 

E 

k=l 



N 



^^(% - Qk-if + eVXwkqk - e^Mk ^ ^fej^j 



(57) 



(58) 



(59) 
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where Wk = 'w{tk) and D^ j = D{tk,tj). The constraints are: qo = xq and = x. We re-write 
GAr(x, t; xo, 0), separating the terms which are constant with respect to the integration variables, 
the hnear terms and the quadratic terms; using a vector notation, we have: 



— 



V 2'nihe 



j exp 



— - [xq + X J + ev awnx — e ax 



/ + 00 
dXexp [-X- AX + 2X- Y] 
-oo 



(60) 



where X and Y are two (A^ — 1) -dimensional vectors defined as follows: 



X 



/ .1 \ 

qN-2 

V QN~1 J 



2he 









V ^ / 



+ 



^ Wl \ 

WN-2 
\ WN-1 ) 
7V-1 



D2,N 
Dn-2,N 



(61) 



We have also used the short hand notation: dX = Jlfc=i '^^^ matrix A is the sum of two 



symmetric, (A^ — l)-dimensional square matrices B and C, whose entries are: 

R - J— n 
he 

Bi,j = , j y^i,i±l, 



im 
We 



Cij — e^XDi j . 



(62) 
(63) 
(64) 

The multiple Gaussian integral of Eq. (j60|) can be immediately evaluated by using the standard 
result: 



dX exp [-X • AX + aX • Y] 

GAr(x, t; xo, 0) then becomes: 

GAr(x,t;xo,0) = 



vr 



N-l 



det(A) 



exp 



-Y • A-^Y 
4 



(65) 



'{m/2ihe)^ 



7rdet(A) 



exp 



Y • A'^Y + 



im 
2he 



,^0 



+ x"^) + e^/Xw]\fX — e'^Xx'^DN^N] ■ (66) 



Note that the integral in Eq. (j65p exists only if A is a positive definite matrix. The result can 
be extended to non-negative matrices, following e.g. the procedure of Theorem 1, page 13 of [3]. 
Accordingly, our results strictly hold only when the correlation function D{t,s) of the noise w{t), 
seen as an integral kernel, is non-negative definite. 

The next step is to take the limit A^ ^ cxd, in which case one encounters two main difficulties: 
the first is to evaluate the inverse matrix A~^, the second is to compute the determinant of A, in 
both cases for any A^. To solve these difficulties, we proceed as in [29!\ . 
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In order to evaluate Y • A~^y, we introduce a twice differentiable function z{s), s £ [0, t], yet to 
be determined, such that, given the vector Z := (zi, . . . zn^i)^ , with Zk = z{tk), k = 1, . . . N — 1, 
the following matrix equation is satisfied: 



AZ = Y. 



For reasons which will be clear soon, we will set z{0) = xq and z{t) = x. We then have: 



(67) 



Y • ^"^Y = Y • Z 



N-l 



im e\/A e^A \ - 

— (XOZI + X ZN-l) + ^ 2^ WjZj - —X 2^ Dj^NZj 



2he 
im 



2^^xo (z(0) + ez'iO) + 0{e')) - —x {z{t) - ez'{t) + 0{e')) + 



WjZj - —X ^ Dj^NZj] 



(68) 



the prime in z'{s) denotes the first derivative of z{s). Inserting the above result in (j66p . we have: 



GAr(x,t;Xo,0) 



{m/2ihe) 
7rdet(^) 



N 



■ exp 



— {xoz'iO) - xz'it)) + '-^Yl ""1^1 + 



N-l 



(69) 



where we have collected all terms of order e or higher, which will vanish in the limit N ^ oo. 
Note that the two terms of (|66p proportional to Xq and x^, which are of order e~^, are canceled by 
the analogous terms in ()68p : this is a consequence of the boundary conditions we have set on the 
function z{s). 

In order to compute the determinant of the matrix A = B-\-C, where B and C have been defined 
in (j62p and (j64p . it is convenient to isolate the determinant of B by writing A = B{I + B^^C), so 
that: 



det{A) = det{B) det(/ + B'^C) 



(70) 



The quantity det{B) will be used to simplify part of the terms under the square root in (|69p . The 



matrix B is equal to B = B with 



Bi^i = 2, 



B 



i,i±l 



■1, Bi^j = 0, J / i,? ± 1. 



(71) 



Let us call the determinant of the matrix obtained from B by removing the first N — k — 1 
rows and columns; the following recursive relation is easy to prove: 



At = 2A 



fc-l — 



A; = 1,... A^- 1, 



(72) 
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where we have set Aq = 1 and A_i = 0. From this relation, by induction one immediately sees 
that Afc = k + 1, which means that det{B) = N. Accordingly we have: 

N-l 



det{B) = n(^ 



m 
2ihe 



(73) 



Collecting all results we can write: 

GN{x,t; xo,0) = 



m 



2i7rhtU]\f{t) 



exp 



— (xo^'(O) - xz'{t)) + ^Y1 ^^-^i + ^(^) 



2h 



, (74) 



where t = Ne, and we have introduced the quantity UAf(t) := det(/ + B ^C). We are now in the 
position to perform the limit N oo; the Green's function becomes: 



G{x,t; xo,0) 



m 



2i7rhtu{t) 



exp 



2^ ,xoz'{0) - xz'it)) + 2 



ds w{s)z{s) 



(75) 



which is defined in terms of the two unknown functions z(s) and u{t) := limTv^oo ""Af (i)- We now 
show how to determine them. 

Written in components and once divided by e, Eq. ()67p becomes: 



im 



N-l ^ 

2^^2 ~ 2zfc + Zk-i) + eA ^ DkjZj = -^Wk - ya^-Dfc.Af, k = I, . . . N - 1; (76) 

i=i 

note that for A; = 1 and for A; = — 1, consistency is assured thanks to the assumption zq = xq 
and zn = x. Taking the limit e — > (A^ oo) we obtain: 



— 2:"(s) + A J dr D{s,r)z{r) = —w{s) 



(77) 



which is an integro-differential equation to be solved together with the conditions z{0) = xq, 
z{t) = X. The above equation determines z{s). In appendix [A] we prove that Eq. ()77p admits a 
unique solution. 

Let us now determine u{t)] we follow once again the procedure outlined in 
Statement 11.11.4, page 451 of [3] one can write 



29[. According to 



det(/ — rjK) = exp 



Tr[i2(^)](i/i 



(78) 



where R{pi) is a matrix satisfying the equation R{ij) = K + fiKR. Applied to our problem, we 
have: 

.0 / 



UN{t) = dei{I + B ^C) = exp 



(79) 
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in our case K = B~^C, thus R{^) satisfies the matrix equation BR(fi) = Ce~^ + fj.CR{fi). Written 
in components, and divided by e, this equation becomes: 

N-l 



im 

2h^ 



[Ri-ij in) - 2Rij (fi) + Ri+ij (fi)] = XDij + /i e ^ A.fc-Rfej (^) , (80) 

fc=i 

with the boundary conditions i?o,j(/u) = Riyj{fj,) = for any j and fj,. Taking the hmit e ^ 0, the 
matrix Rij(fj,) becomes a function i?(s,r, /i) satisfying the following integro-differential equation: 



q2 ft 

R{s,r,fi) - Xfi / dlD{s,l)R{l,r,n) = XD{s,r), (81) 
Jo 



which has to be solved together with the boundary conditions R{0, r, /i) = R{t, r, /i) =0 for any 
r G [0,t] and /x G [—1,0]. The existence and uniqueness theorem of appendix[X] applies also to this 
equation. In the limit N — > oo, Eq. (j79p becomes: 

rO / rt 



u{t) = exp 



J dsR{s,s,fi) 



(82) 



Eq. ([75]) . with z{s) defined by ([77|) and u{t) defined by ([82]) . gives a complete description of the 
Green's function of Eq. (|53|) . This is the main result of this section. 

The function z{s) solution of Eq. ([77]) depends on the end points xq and x because of the 
boundary conditions. This is not evident in the expression (j75p for the propagator. In order to 
make such a dependence on xq and x explicit, we rewrite z{s) as follows: 

z{s) = f{s)xo + g{s)x + h{s), (83) 

where /(s) satisfies the homogenous integro-differential equation: 

'-^f\s) + xj\rDis,r)fir) = 0, (84) 

with boundary conditions /(O) = 1, f{t) = 0, while g{s) solves the same equation but with 
boundary conditions g{0) = 0, g{t) = 1. The function h(s) instead satisfies the non-homogenous 
equation: 

^ h"{s) + ^l^dr D{s, T)h{r) = ^w{s), (85) 

with boundary conditions /i(0) = h{t) = 0. Note that /(s), g{s) and h{s) depend also in the 
parameter t, so one should more properly write ftis), gt{s), ht{s); sometimes we will omitt the 
pedex t when no confusion arises. Here and in the following, the prime denotes differentiation with 
respect to the variable in parentheses. 
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With the help of (j83p one can write the Green's function of Eq. ()75p in the following way: 



G{x,t; xo,0) = ■\J 2i7rhtu{t) ^"^^ -Atxl - Atx'^ + BtXQX + CtXo + Vtx + £t , (86) 

where the dependence on xq and x has been made explicit. The coefficients At, At and Bt are 
deterministic functions of time, defined as follows: 

At = kfliO), At = -kg'tit), Bt = kinit)-g'tiO)), k = —, (87) 



2h 



while the coefficients Ct, Vt and £t are random processes: 



Ct = -kK{0) + ^ [ w{l)ft{l)dl, (88) 

Vt = kh[{t) + ^ f w{l)gt{l)dl (89) 
^ JO 



= V / w{l)ht{l)dl. (90) 

^ JO 

The propagator (j86p takes a simpler and more symmetric form when the noise is time translation 
invariant: D{t,s) = D{t — s); according to the symmetry in the (t, s) variables one also has: 
D{t — s) = D{\t — s\). In this case, one can easily verify that if f{s) is solution of Eq. (1840 with 
the boundary conditions /(O) = 1, f{t) = 0, then f{s) = f{t — s) is still a solution Eq. ([84l) . with 
boundary conditions /(O) = and f{t) = 1; in other words: f{s) = g{s). Then (j83p becomes: 

z{s) = f{s)x^ + f{t-s)x-h[s), (91) 

and the Green's function (j86p simplifies as follows: 



G(x, t; xo, 0) = ^j ^i-Khluit) "^""^ ["'^^(^o + a;^) + -Si^oa; + CtXo + Aa; + Et] , (92) 
where the deterministic coefficients now are: 

At = kfliO), Bt = 2km, (93) 

and the stochastic terms read: 

Ct = -kh'M + ^ f dlw{l)ft{l), (94) 

^ JO 



A = kh't{t) + ^l^ 



dlw{l)ft{t-l), (95) 



ft = ^ f dlw{l)ht{l). (96) 

^ JO 
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Accordingly, for a general time translation invariant noise, the structure of the Green's function is 
the same as in the familiar white-noise case 



101, 



30l |. as we will review in Sec. El Note that due 



to linearity, the general solution of Eq. (I85p can be written as h{s) = h^{s) + h^{s), where h^{s) 
solves the homogeneous equation and h^{s) is a particular solution of Eq. (j85|) . Often the 

particular solution admits an integral expression of the form: 



dl r{s, I), 



(97) 



where r{s,l) is some function. Eq. (j97p implies that h^{0) = 0. One can immediately verify that, 
by taking h^{s) = —h^{t)f{t — s), the function: 



h{s) = h^{s) - h^{t)f{t - s) 



(98) 



solves Eq. (jSSp . with the correct boundary conditions. Thus also h{s) can be expressed in terms 
of /(s), once hF{s) is given as in ([g?]) . 

A concluding remark is at order. In this paper we have computed the Green's function 
G(x,t;xo,0) by exploiting the path-integral formalism. The same result can be obtained also 
by resorting to the standard operator formalism, starting with the expression for the propagator 
first proposed in 



25|. 



IV. THE NON MARKOVIAN STOCHASTIC SCHRODINGER EQUATION 



The result of the previous section allows us to rewrite Eq. (j53p in a simpler form, removing 
the functional derivative which is not easy to handle. In this way, we will prove the correctness of 
the Ansatz introduced in Sec. 4.3 of Ref. [l7], according to which the functional derivative can be 
replaced with a suitable function of the position and momentum operators. 

The starting point is Eq. (j56p where the solution (/>t(x) of Eq. (j53p is expressed in terms of 
the Green's function G(x,t;xo,0). Using standard path- integral techniques, e.g. generalizing the 



procedure described at page 509 of 



3ll |. one can show that G(x,t;xo,0) given by (f54|l . with S[q\ 



given by (j55]) . satisfies the following equation: 



d_ 
di 



G{x,t; xo,0) 



I p 
h 2m 



G{x, t; xo, 0) + G{x, t; xq, 0) , 



(99) 



with: 



G{x,t; xo,0) :-- 



l-q{t)=x i-t 

-2Xx / V[q] / dsq{s)D {1,8)6^^"^ 

Jq{0)=xo Jo 



(100) 



18 



The path integral in (jlOOp cannot be trivially reduced to because q{s) appears, which de- 
pends on the entire time interval [0,t]. Thus we have to re-calculate G{x,t; xq,0) as we did with 
G{x, t; xq, 0). Proceeding as in Section HIH we define G{x, t; xq, 0) = liniTv^oo Gn{x, t; xq, 0), where 
Gn{x, t; Xq, 0) is: 



GNix,t;xo,0) 



-2Xq 



m 



exp 



im 
2he 



Cq -|- x^) -|- eV\w!s[X — e^Xx^D 



N,N 



r+oc ( \ 

/ dX e^KTx + 0(e) exp[-X-AX + 2X-Y] 



(101) 



where 0(e) collects all terms of order e, = (0, . . . -Data,., ... 0) where D^^ is placed in the fc-th 
position, and X, Y, A are defined in Eqs. (I6ip and (|62p - (|64p . Applying the following generalization 
of Eq. dMD: 



dX K'^X exp [-X • ^X aX • Y] = 

and taking the limit for e ^ one finds that 



"-K"^A-iY exp 



det(^) 2 



Y • A~^Y 



(102) 



^(x, t; xo, 0) = — 2Ax / 2;(s)D(t, s)G(x, t; Xq, 0) , 



(103) 



where z(^s) is defined in Eq. (I77p . Thus, using (I99p and ()103p one can formally write (x = q, the 
position operator): 



di 



G{x,t;xo,0) 



I p 



+ \^qw{t) - 2\q / dsz{s)D{t 



G(x,t;xo,0) . 



(104) 



h 2m 

Note that it is possible to find this same result by inserting G(x,i;xo,0) directly into Eq. (j53p . 
However, to reach the result of Eq. (|lU4p . one encounters a term of the form 5z' [l) / 5w{s) which 
can be simplified only under the assumption that the time derivative and the functional derivative 
acting on z{s) commute. This is not true in general, and one has to check it by finding the explicit 
functional dependence of z{s) on w{s). This fact makes this second approach more cumbersome 
than the one we have followed. 

Eq. (|104p is not written in a useful way yet, since z[s) depends both on xq and x. We make 
explicit this dependence on xq and x = qhy resorting to ([83]) : 

ft 



G(x, t; Xq, 0) 



-2Ax 



ds D{t, s) (/(s)xo + g{s)x + h{s)) 



G(x,t;xo,0) . 



(105) 



The term in (jlOSp depending on xq does not allow to perform the integration over the xq variable, 
to compute (pt from (j56p . This term can be quite easily written in terms of the position operator 
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q and the momentum operator p = —ih-^ acting on G{x, t; xq, 0), as follows: 



xoG{x,t;xo,0) 



9'{o) - m 



2ih d ^, 

— — G(x,t;xo,0) 
m ax 



+ l2xg'{t) + h'{t) + ^-^^j^ dlwil)g{l) \ G(x,t;xo,0) 



Collecting all the above results, we then find that (j)t{x) satisfies the following equation: 

Mx) , (106) 



I ^ + ^qw{t) - 2Xq [ dsD{t,s){qat{s) +pbt{s) +ct{s)) 



where 



2 m 



bt(s) 



ct{s) = h{s) + — 



mg'{0)-f'{t) ' 
fis) 



(107) 
(108) 

(109) 



This is the non-Markovian stochastic Schrodinger equation, whose Green's function is given by 
Eq. (j86[) . When the noise is time-translation invariant {D{t,s) = D{\t — s\)), the functions at{s), 
his), ct{s) simplify as follows: 



f (0) 

atis) = f{t-s) + Ly.f(s), 



his) 



ct[s) 



1 fjs) 
m fit) ' 



(110) 
(111) 

(112) 



As anticipated, Eq. (jl06p agrees with the ansatz first proposed in [ir]. 



V. WHITE NOISE LIMIT OF THE GREEN'S FUNCTION 



Having found the Green's function associated to Eq. (|53p . the first task is to check that, in the 
white noise limit, it coincides with the propagator already known in the literature [lol . 3o| . Since 
D{s,r) = 6{s — r) is time translation invariant, we can use the expression (j9ip for z{s). In the 
white noise case, Eq. ([8l|) reduces to the second order differential equation: 



/"(«)- ^/(s) = 0. 



m 



/(0) = 1, f{t) = ^, 



(113) 
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whose solution can be easily computed: 



/(s) = coshfs — cothut sinhus , 

l + i ^ [hX 

V := -^uj, uj := 2W — 
2 V m 



(114) 
(115) 



Eq. ([85]) for h{s) instead reduces to: 

2iXh 



h"{s) 



ih 



m 



h{s) = \/Ail'(s), 



m 



h{0) = h{t) =0. 



(116) 



Using once again time translation invariance, we know that h{s) takes the form ()98p . where a 



particular solution h^{s) is 



h^(s) = \/A / dr wlr) smhv(s — r). 

mv Jq 



(117) 



The Green's function then has the form (j92p . where the deterministic coefficients At and Bt are 
equal to: 



At 



A 



coth vt . 



2A , _i 
Bt = — smh vt, 

V 



(118) 



while the random coefficients Ct, T>t and £t are equal to: 



sinh vt 
^ sinh vl 
sinh vt ' 



rr „ ,,.smhv(t — l) 
Ct = VX dlw{l) ^ ' 

Jo 

Vt = VX [ dlw{l) 
Jo 

£t = —( f dl w(l) sinhvl [ dsw(s)— ^ ^ 

4 Jo smhvt 

dlw{l) I dsw{s)smh.v{l — s)\ . 
Jo J 

After some manipulation and an integration by parts, one can rewrite £t as follows: 

^ /"* „ sinhu(t - /) , , 

Ct = — dlwil) ; / dsw{s)s\mivs 

2 Jo ^ ^ smhvt Jo ^ ' 



(119) 
(120) 



(121) 



2 rt \^ dsw{s) sinh vs 



V 

^ JO 



sinh^ vl 



(122) 

(123) 
(124) 



When £t is written as in Eq. (jl22p . then the Green's function (j75p reduces to that of Eq. 6.2.42, 



page 180 of 



30l |. When £t is rewritten as in Eq. ()123p . it coincides with the corresponding term 



of 
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We now compute R{s,r,iJ,), which in turn determines u{t). In the white noise hmit, Eq. (j81|] 
becomes: 

a2 



R{s,r,fj,) + v fiR{s,r,fj,) = —v 5{s — r), 



(125) 



with R{0, r, /i) = R(t, r, fj,) =0 for any r, ^. R{s, r, //) clearly displays a discontinuity on its first 
derivative for s = r. The solution is found by imposing on the solution of the homogeneous equation 
the correct discontinuity conditions on the first derivative, together with the boundary conditions. 
R{s,r,fj,), having the correct boundary conditions in s = and s = t, is easily determined: 



a sinh. iv^/J2s 



s < r . 



Ris,r,n) = < 

b sinh iv ^/JI{s — t) s>r, 
where a and b are two constants determined by the following relations: 

a sink iv^/JIr = bsm]iiv^/J2{r — t) , 



(126) 



(127) 



2i 

mv 



b^/Ji cosh iv^/Jl{r — t) — a^/Jl cosh iv^/Jlr 

the first equation represents the continuity of R{s, r, /x) in r and the second equation encodes the 

discontinuity of the first derivative. 

After some algebra on finds the following result: 

' iv sinh iv y/JI{r — t) sinhiv^^s 



Ris,r,n) = < 



sinh ivyfjit 
iv sinhivyfflr sinh iv yfji^s — t) 



sinh iv^/Jlt 



s < r, 



s > r. 



(128) 



One then has: 



dfi / dr R{r, r, jj) 



In. 



siuhiv^JJlt 



(129) 



the above formula, when evaluated in /i = — 1, produces an extra i factor which cancels the one in 
the argument of sinh. Substituting this result into Eq. ([82|) . one finds that: 

sinh vt 



u{t) 



vt 



Substituting Eq. (jl30p into (j75p the prefactor in front of the Green's function becomes: 



A 



■nv sinh vt 



which is the prefactor known in literature 



lOl, 



3d, 



(130) 



(131) 



33] 
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VI. EXPONENTIAL CORRELATION FUNCTION 

The results of section IIIII are important by themselves. However, they acquire a further impor- 
tance when Eqs. (j77p and (I8ip can be solved for some specific, physically meaningful, correlation 
function. In this section we treat the case in which the correlation function has an exponential 
damping, i.e. 

D{t,s) = le-^\'-'K (132) 
where 7 is the inverse of the correlation time. This type of correlation function was first considered 

n 

in [171]. Since D{t,s) is time translation invariant, we can use the expression (|9ip for z{s). The 
equation for f{s) reads: 

_n^j\r e-^l«-H/(r) = , (133) 
with oj defined (I114p . The equation has to be solved with the boundary conditions /(O) = 1, f{t) = 



34( 1 : one differentiates 



0. The solution can be found using the following mathematical procedure 
the equation twice, obtaining: 

/""(.) + ^^Wm + 7^ (-^ 1^ dl e-^(-0/(;) -'^1^1 e^(-')/(/)) = ; (134) 

next, using once again Eq. (I133p . one replaces the third term with — 7^/"(s). In this way Eq. (jl33p 
is transformed into the following fourth-order differential equation: 

/""(,) _^Y'(s) + i7V/(s) = 0. (135) 

The general solution of this equation, which has the advantage of not having an integral term, is 
easily found and can be expressed as a sum of hyperbolic sines and cosines as follows: 

f{s) = fi sinhi;is -|- /2 sinhu2S + /s cosh vis + fi cosht;2S , (136) 

where fi are constants, to be determined by the boundary conditions, while Vi are the two non- 
symmetric roots of the bi-quadratic characteristic polynomial associated to Eq. (jl35p . i.e. 



vi,2 = y\{i^±C), C = Vi^- 4i7'^' ; (137) 

the other two solutions are: ^3 = —vi and ^4 = —V2- It is easy to verify that in the white noise 
limit (7 ^ cxd): — > 0(7) and V2 ^ v defined in (jll4p . 
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Since Eq. (jl35p is a fourth-order differential equation, to determine a specific solution we need 
four boundary conditions. Two of them are /(O) = 1 and f{t) = 0. In order to find the other two 
conditions, we proceed as in [34] • We evaluate Eq. (1133P in and in t: 



no) 



fit) 



2 rt 



dle-^^f{l) = 



(138) 



2 rt 



and we replace f{l) in each expression, using Eq. (jl35p : 



f"{0) + - / dle'"''if""{l)-^'f"{l)) = 
^7 Jo 



f"{t) + ^ [ dle~^('~^Hr{l)-j'f"{l)) = 0. 
27 Jo 



(139) 



Integrating by parts these two equations, after some algebra one finds that the complete set of 
boundary conditions is: 

'/(o) = i, r(o) = 7r(o) 

(140) 

^ /(t) = o, r{t) = -jf"it). 

The solution of this system involves a long calculation and leads to a complicated expression 
for the coefficients fi in (jl36p . In order to write the solution /(s) in a compact way, we introduce 
the following new coefficients: 



3„,3 



f 1 W 



dk 



l'^2 ) 



(141) 
(142) 
(143) 
(144) 



where k = 1,2 and A; = 2 if A; = 1, A; = 1 if /c = 2. With the help of these prescriptions, /(s) can 
be written as follows: 



/(s) = Ms) 



J2k [''t smhvk{t — s) + coshffc(t — s) 



where: 



J2k [2c + sinh Vkt + cosh Vkt] 



az cosh vzt + bz sinh vzt 



dk sinh vj^t — c cosh vj^t . 



(145) 

(146) 
(147) 
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Resorting once again to time translation invariance, we know that h{s) takes the form 
when the particular solution h^{s) has the form (I97p . In order to find h^{s), we have to consider 
the non-homogeneous equation: 

h"(s)-'-^ t dre-'\'-Ms) = --^w{s). (148) 
2 Jo 

Following the same calculations which brought from Eq. (jl33p to Eq. (jl35p . one finds the following 
fourth-order non-homogeneous equation for h{s): 



h""{s) - ^\"{s) + ij^oj^h{s) = iw"{s) - ^Ms)) ■ (149) 



A particular solution of this equation, having the form (|97l) . is 



h^(^s) = / /,(/) {w"il) - jMI)) dl , (150) 



m 







where /s(/) is the solution of Eq. (jlSSp . thus having the general form (jl36p . with boundary con- 
ditions: /s(s) = f's{s) = f's'is) = 0, f's"{s) = 1. After some simple calculations, one finds the 
following expression: 

We now have all functions which are necessary to compute the Green's function. This takes the 
form (j92p . where the coefficients At^St are defined in Eqs. (j93p - ()96p in terms of the following 
functions: 

jf^Q^ ^ Efc ^fc {4 cosh Vkt + u'l sinh Vkt + 4) ^^^^^ 
* J2k [2c + sinh Vkt + v!l cosh Vkt\ ' 

^/^^^ ^ Efc [rt + 4 cosh Vkt - c sinh Vkt) ^^^^^ 
J2k [2c + sinh -|- cosh f^t] ' 

and: 

/i^'(O) = (154) 
^p,^^^ ^ ^^/A^ /■*^(_i)fccoshi;fe(t-0 K'(0 -7^w^(0) • (155) 

In terms of these functions, the entire dynamics of the free particle can be analyzed. 

Eq. ()8ip for R{s,r,fi) can be solved similar to Eq. ()77p . One then finds the explicit expression 
of the prefactor of the Green's function (|92l) . This goes outside the scope of the present paper, so 
we omit the calculation. 
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VII. DENSITY MATRIX EVOLUTION AND IMAGINARY NOISE TRICK 



One of the major difficulties connected to Eq. ([6]) arises when one wants to compute average 
values of the form Ep[(V't |0|^t)] = E(q[((/)( |0|i;^)t)], where (pt solves Eq. ([6]) and ipt is the corresponding 
normalized solution; O is a generic self-adjoint operator. Averages of this kind are particularly 
important, as they represent physical quantities, directly connected to experimental outcomes. The 
difficulty in computing such averages lies both in the difficulty in solving Eq. ([6]) and in the fact 
that (pt depends non-trivially on the noise w{t). 
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3g] allows 



In the white noise-case, a very helpful trick, known as the imaginary noise trick 
to simplify considerably the problem. Let us consider the following class of stochastic differential 
equations: 



N 



N 



n=l 



n=l 



4^ 

(156) 



where C = + i^i is a constant complex factor. By using standard Ito calculus one can show that 
Ep[(V'f lOI^/'f )] depends only on |^|^: there is, so to speak, an invariance under phase change in the 
coupling constant. When = 1 one recovers Eq. ([1]), while when = i, one obtains a standard 
Schrodinger equation with a stochastic potential, generating a unitary (thus non-collapsing) evolu- 
tion, without non-linear terms; this equation is much simpler to analyze than Eq. ([T]), nevertheless, 
at the statistical level, it gives the same results as Eq. ([1]) does. 

In [2l| it has been shown that the imaginary noise trick holds also in the non-white noise case, 
at least to second order in the perturbation expansion with respect to the parameter \/A. Aim of 
this section is to prove that this is an exact property which holds to all orders. 

Let us consider the following class of equations: 



-'-H + iV\qw{t) - 2^^VXq [ dsD{t,s)-^^ 
n Jq dw{s) 



(157) 



obviously, when ^ = 1 one recovers Eq. for a free particle. Since Ep[('i/'t |0|?/'t)] = 

EQ[{(j)t\0\(j)t)] = Tr[Opt] where the density matrix pt is defined as: pt := Ep[|'i/'t)(V'i|] = IEQ[|<At)('At|]) 
in order to prove the required property it is sufficient to show that pt depends only on The 
propagator associated to Eq. ()157p reads 

rq(t)=x 



G^{x,t; xo,0) 



q{0)=xo 



where 



sdq] 



ds 



-So[q] + CV\qis)w{s) - C> 



Jo 



dr q(r)D(s, r) 



(158) 



(159) 
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and So[q] is the standard action associated to the quantum Hamiltonian H. The propagator 
J(x, x', t; xq, Xq, 0) associated to pt{x,x') := {x\pt\x') can then be expressed as fohows f28i ]: 



J{x,x',t;xo,XQ,0) = Eq [G^(x,t;xo,0)G|(x',t;xo,0)] , 



(160) 



where * denotes complex conjugation. Substituting the definition ()158p for the wave function's 
propagator into (jl60p . and exchanging the path-integration with the stochastic average, the stochas- 
tic terms average as fohows: 

rt 



exp 



exp 



1 



t rt 



JO 



ds / drD{s, r) {eHs)qir) + e^Wis)q'{r) + Xq{s)q' (r) + \^nq{r)q' (s)) 



(161) 



The propagator then becomes: 
J(x,x',t;xo,Xo,0) = 



q(t)=x rq'(t)=x' 

V[q] / Pte']exp 

g{0)=xo Jq'{0)=x[, 



d^{ji^o[q]-'j^So[q'] 



f drD{s,r) {q{s) - q'{s)) [q{r) - q'{r)) 



(162) 



which depends only on We can conclude that the evolution for the density matrix is indepen- 
dent from a phase change in the coupling with the noise. 



A. Average values of physical quantities 



As anticipated, this phase change invariance provide a very handy tool to compute average 
values, such as the mean position and mean momentum. In place of Eq. (j53p . let us consider the 
equation: 



I p . I— , , 
2m ^ ^vAg'(y(t) 



0i(x) 



(163) 



which belongs to the class ()157p . with £, = i- This is a linear, unitary, norm-preserving standard 
Schrodinger equation for a free particle under the influence of a stochastic potential. The evolution 
of the stochastic average of the mean value of an observable O is given by the following equation: 



m{[P^O])t]-t^EQ[w{t){[q,0])t 



(164) 



where, as usual, {•)t = {4>t 
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The mean in position and momentum can now be easily computed. Remembering that {(j)t\4>t) = 
1 and that the mean of w{t) is 0, one finds: 

^Eq [{q)t] = Mt] , (165) 

^EQ[(p)t] = 0; (166) 

as we see, we recover Newton's equations for a free particle. In particular, also in the non-Markovian 
case, like in the white-noise case, the momentum of an isolated system is conserved, in the average. 
In sec. I Villi where we analyze the time evolution of Gaussian states, we will see that the fluctuations 
around the average are inversely proportional to the mass of the system. These two facts, together 
with the fact that the collapse scales with the size of the system, lead to the following result: the 
wave function of an isolated macro- object behaves, for all practical purposes, like a particle moving 
deterministically in space according to Newton's laws. 

Because of the many experimental implications [3], it is important to check how the mean free 
energy Hq = p^ /2m evolves in time. This is most easily computed by shifting to the Heisenberg 
picture, in which case one finds 



2lC: 

-EQ[(Fo)t] = — / dsD{t,s)- (167) 
at m Jq 



this is the expected generalization of the well-known white noise formula for the mean energy 
increase in collapse models [q. [l^. [l^. For a physically reasonable correlation function such as the 
exponential one of Eq. (I132p . one has: 

A/i^ / R-^* - 1 \ 
Eq [{Ho)t] = {Ho)o + ^(t + — ■■ (168) 

as we see, also in this case the mean energy increases linearly in time, without reaching a steady 
value. More generally, let us assume time translation invariance, and let us consider the spectral 
decomposition of a generic correlator D(\t — s\): 

/•oo 

D{\t-s\) = / dujD{uj) COS u; {t - s); (169) 
Jo 



we have: 

rt /•oo 
'0 

The above formula shows that the energy production is nonzero also at large times, unless the 
correlator has a cutoff at low frequencies [2l[ |. 



/' 

Jo 



dsD{t,s) = / —D{u/t)smu > tt-D(O). (170) 

Jq U t^+OO 2 
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An important lesson to learn from the above analysis is that non-Markovian terms do not intro- 
duce thermalization effects in the evolution of a quantum system; such effects can be introduced 
only by modifying the form of the operator coupled to the noise, as first discussed in jsj]. Never- 
theless, non-Markovian terms are extremely important, as they affect the time evolution of physical 
quantities, thus t he p redictions of collapse models, at small times. A significative example has been 
first provided in 3a]. 



VIII. EVOLUTION OF A GAUSSIAN WAVE FUNCTION 



In this section we study the time evolution of Gaussian wave functions, whose form is clearly 
preserved by the Green's function (j86p : for simplicity, we will assume time translation invariance, 
so that ((86]) reduces to (f92]) . A generic Gaussian state of the form: 

Mxo) = 

evolves in time to the following state: 

Mx) = 

where: 

at = 
Pt = 
It = 

and the functions At - £t have been defined in Eqs. (f93l) - (i96]) . We can draw some important 
conclusions about the time evolution of Gaussian wave functions: 

• Since neither At nor Bt depend on the noise w{t), the time evolution of the spread of the wave 
function is deterministic, as in the white noise case Q]. An initially spread-out Gaussian wave 
function shrinks in space, reaching an asymptotic final spread. In the case of an exponential 
correlation function, we can provide the explicit expression for the asymptotic value for the 
spread (see next subsection). 

• By making the substitution: 

X , (176) 



exp[-aoa:;o + /^oa^o + 7o] 



(171) 



exp[-afX^ + Ptx + 7t] , 



(172) 



At 



A{ao+At) 
Ct + /?o 



^, , . .Bt + Vt. 
{Ct + Pof 



70 + <St + 



4(ao + At) ' 



(173) 
(174) 
(175) 
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one can easily prove that Eq. ([53|) . thus the propagator ([92]) . does not depend on the mass 
ra of the particle. Thus, one way to see the effect of the mass on the global dynamics is to 
take the evolution of the wave function for the reference mass mo, and then "shrink" the 
space coordinates by a factor y^^^o/^- This leads to the amplification mechanism, which 
is the characteristic feature of collapse models: the bigger the system, the faster the wave 
function shrinks in space. 

• The center {q)t = /2a^ of the Gaussian wave function (jl72p evolves randomly in space, 
as expected. Its average value has already been computed in (jl65|) . and evolves classically. 
Moreover, due to the independence of the dynamics from m under the substitution (jl76p . 
one can conclude that the fluctuations of {q)t around the classical motion go like m~^/^ since 
the following equality for the variance Vg := \J^q [{<l)t — EQ)[(<?)t]]^ follows immediately 
from (fTTHIl : 




(177) 



where the apex m indicates with respect to which mass the variance is computed. Eq. (|177p 
in agreement with the white noise case. This means that the bigger the system, the smaller 
the fluctuations: in the microscopic case, in this way, one recovers classical determinism for 
all pratical pourposes. 

• The mean momentum {p)t = h[(3l — {a\/ a^)0f] also evolves randomly in space. Its average 
value is constant in time (see Eq. (|166p l. as expected from a free particle, while the fluctu- 
ations around the average increase like m}!'^ . If however we consider the fluctuations of the 
mean velocity, we have that they decrease like mr^l'^ . Thus also in this respect one recovers 
classical determinism at the macroscopic level. 

The above remarks show that the non-Markovian QMUPL model shares all the important features 
of the white noise model. More quantitative details can be given by taking a specific expression 
for the correlation function D{t, s), as we will do in the next subsection. 

A. Exponential correlation function 



We study in some more detail the time evolution of the Gaussian solution (|172p . when the 
correlation function is exponential as in (jl32p . We focus our attention on the spread in position 
a{t) which is given by the inverse of twice the square root of the real part of Qt, whose analytic 
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tr(f) 
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2x10"-* 4x10"" 6x10"^* SxlO"'-* 1x10"^^ 



FIG. 1: Linear/linear graph showing the time evolution of (T{t) for a 1 Kg particle, for different values of 
7, for the same initial condition (t(0) — 1 and with Aq = 5.00 x 10^'^ m~^ sec'^. The curve with 7 — 00 
corresponds to the white noise case. Time is measured in seconds, distances in meters. 




100 1000 10* 10' lo' 10' 



FIG. 2: Log/linear graph, showing the large time behavior of <j{t) for a 1 Kg particle, for different values of 
7, for the same initial condition <j{Q) = 1 and with Aq = 5.00 x 10^"^ m^^ sec^^. The white noise case would 
appear in the graph as a straight line with value 1.27 xlO~^^ m. Time is measured in seconds, distances in 
meters. 

expression is given explicitly by Eqs. (|173p . ([93j) and ()152p . Figs. [1] and [2] show some typical 
behavior of a{t) for small times (Fig. [T]) and large times (Fig. [2]), for different values of 7 and for 
m = 1 Kg. As we see and as we expect, the larger the value of 7, the stronger the noise and the 
faster the collapse in space. 

Fig. [3] shows how the spread of a Gaussian wave function associated to a particle having a mass 
m = 1.01 X 10^'^ Kg (which is the total mass of a system of 6.02 x 10^'^ nucleons), having initial 
spread o"(0) = 1 m, decreases after 10^^ sec, as a function of 7. The value of the coupling constant 
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o-(t=10-^) 




7 



FIG. 3: Linear/linear graph, showing to which value an initial spread cr(0) = 1 m decreases after 10~^ sec, 
as a function of 7. The mass of the particle has been set equal to m = 1.01 x 10^'^ Kg, corresponding to the 
total mass of a system containing an Avogadro's number of nuclcons. The coupling constant Aq has been 
given the GRW- value 5.00 x lO^'^ m^^ sec^^. Larger values of 7 imply a faster collapse in space. Time is 
measured in seconds, distances in meters. 



Aq is the GRW value ([3]) . GRW have chosen the value of the coupling constant in order to ensure 
that the wave function of a system of at least an Avogadro's number of particle collapses, within 
a time interval of 10"'^ sec, below 10""^ m. Fig. [3] shows that also for relatively small values of 7, 
the non-Markovian collapse models preserves this same feature of the white-noise model. 

jig. [3] displays the same graph as that of Fig. [31 with the typical values chosen by Adler. 



In [l3|] Adler sets the value of the CSL-coupling constant 7csl equal to 2 x 10^^^^^ cm^ sec~^ 
by noting that in the process of latent image formation, which has a characteristic time of about 
3.33 X 10~^^^ sec, a number of atoms approximately equal to 20 is involved. He then assumes that 
the collapse process ensures that the reduction occurs already at this stage, from which the value 
7csL = 2 X 10^^^^^ cm^ sec^^ comes. The graph shows that even for relatively small values of 7 
(Adler and Ramazanoglu [s^ have shown that choosing 7 ~ 10^^ sec~^ already changes significantly 
the predictions of non-Markovian models with respect to white-noise models), an initially spread- 
out wave function shrinks rapidly below 10"'' m, which corresponds to a well-localized wave packet. 

For an exponential correlation function, we can write down the analytic expression of the asymp- 
totic value of at, which is: 

ITTl 

Ooo = lim at = -7^(f 1 + V2-7), (178) 
t^co zn 

with vi^2 defined in (|137p . From this, the asymptotic spread in position can be easily computed. 
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FIG. 4: Linear/linear graph, showing to which value an initial spread (t(0) = 1 m decreases after 3.33 x 10^^ 
sec, as a function of 7. The mass of the particle has been set equal to m = 1.06 x 10~^* Kg, which is the mass 
of a system of n^N (n — 5640, N = 20) nuclcons, with n and N taken from Eq. (8) of [l3(. The coupling 
constant Aq has been given Adler's middle value 1.12 x 10^ m^^ sec~^. The dashed line corresponds to the 
white noise case. Time is measured in seconds, distances in meters. 



In particular, if we take the white noise hmit 7 ^ cxd, one obtains the finite value: 



which matches with the value known in the hterature 



a. 



IX. CONCLUSIONS 



We have thoroughly investigated the dynamics of a free quantum particle as described by the 
non-Markovian QMUPL model of spontaneous wave function collapse. We have provided an explicit 
formula for the Green's function; we have shown that it reproduces the well-known white noise case, 
and have analyzed the physically interesting case of an exponential correlation function. We have 
proven that the non-Markovian model shares all the most important features of the corresponding 
white-noise model; we have described in particular the evolution of Gaussian wave functions. 

There are of course several other important issues which need to be investigated. In particular: 

• It is important to set on a rigorous mathematical ground the change of measure defined 
in Q, and to derive the analogous of Girsanov's theorem which holds for the white- noise 



case. A sketch of the proof can be found in Appendix B of Ref. 



3- 



It is also important to set a bound on the spread of the general solution as a function of time, 
in order to see how it decays in time. Such a formula would be crucial for setting a lower 
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bound on the value of A by guaranteeing that the model collapses sufficiently big systems in 
sufficiently short time. 

• In the case of an exponential correlation function, it would be interesting to prove if any initial 
state collapses asymptotically to a Gaussian state whose spread in position and momentum 
is fixed and given by (jl78p . In this paper we have proved that this is true only for the special 
case of initially Gaussian states. A similar general theorem has been recently proved for the 
white-noise case {lo| . 

These problems will be the subject of future investigation. 
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finally wish to thank A. Fonda for his guidance through the proof of the existence and uniqueness 
theorem of Appendix lAl 

APPENDIX A: AN EXISTENCE AND UNIQUENESS THEOREM 

In this appendix we prove that equation (i77|) : 



with boundary conditions z{0) = xq, z(t) = x, admits a unique solution. The same theorem applies 
also to Eq. (jSip for R{s,r,fj,). In order to simplify the proof, it is convenient to set both boundary 
conditions to zero. This can be done without loss of generality, as follows. Let us define the new 
function: 



Acknowledgments 




(Al) 



z{s) = z{s) - y{s) 



X — Xq 



s + xo; 



(A2) 



t 



obviously, z(0) = z{t) = 0. Moreover, z{s) solves the following equation: 




(A3) 



where 




(A4) 
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It is then sufficient to prove the following theorem. 

Theorem. Let D{t, s) be a real continuous function on [0, t] x [0, t], symmetric in its two variables; 
let /(s) belong to C[0, i]. Then the integro-differential equation 

iz"{s)+ [ drD{s,r)z{r) = f{s), (A5) 
Jo 

with boundary conditions 2:(0) = z{t) = 0, admits a unique solution z{s) £ C^[0,t]. 

Proof. Let C C'^[0,t] the set of functions z{s) G C^[0,t] such that z{0) = z{t) = 0. Define 

the following operators: 

D : Pp) ^ C7[0,t] , D[z{s)\ = iz" {s) , (A6) 

I : C[0,t] ^ C[0,t], I[z{s)] = [ drD{s,r)z{r); (A7) 

Jo 

with these definitions Eq. ()A5P becomes 

D[z{s)]+I[z{s)] = f{s). (A8) 



The integral operator I is compact (see, e.g., Theorem 8.7-5, page 454 of 
and its inverse reads: 



39(1). while D is invertible 



D-^ : C[0,t]^V{D), (A9) 
= -i du £ dv g{v) - \ du £ dv g{v)j . (AlO) 



Since it has an integral form, also D ^ is compact. We can write Eq. (jASP as follows: 

(Id - D-H)[zis)] = f{s) , z(0) = z{t) = , (All) 

where Id is the identity operator and /(s) = D^^f{s) G C^[0,t]. 

Our problem reduces to showing existence and uniqueness of solutions for Eq. (jAlip . Since the 
product of two compact operators is compact, then also D^^I is compact. We now use Fredholm's 



Alternative Theorem (Theorem 8.7-2, page 452 of 39|]), according to which in order to prove the 



theorem it suffices to show that the homogeneous equation associated to Eq. (jAlip . i.e. 

(Id - D-^I) [z{s)] = , z(0) = z{t) = , (A12) 

admits only the trivial solution z{s) = 0. 

To prove this, we separate the real and imaginary part of z{s) which we respectively denote 
with the superscripts r,i, obtaining 

(z")"(s)+ / drD{s,r)z\r) = 0, (A13) 
Jo 

{z^)"{s)- [ drD{s,r)z'^{r) = 0, (A14) 
Jo 
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with boundary conditions z^{0) = z^'{t) = z^{0) = z^{t) = 0. Multiplying these two equations 
respectively by z^{s) and z^{s), and integrating by parts one finds: 

ds [(z^)'(s)]^ - fdsz^'is) [ drD{s,r)z\r) = 0, (A15) 
Jo Jo 

[ ds [(z')'(s)]^ + / dsz'is) [ dr D{s,r)z''{r) = 0. (A16) 
Jo Jo Jo 

Exploiting the symmetry of D{s, r) in its variables, we can sum the two equations, obtaining: 

' ds[{{z-y{s)f + {{z')'{s)f]=G; (A17) 



this implies that z^'{s) and z^{s) are constants and, applying the boundary conditions, these con- 
stants are equal to zero. We can then state that Eq. ()A12p admits only the trivial solution z{s) = 0. 
This proves existence and uniqueness of solutions for that Eq. (|Alip . and thus for Eq. (f77I) . 
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